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Simultaneous observations of the quasar 3C 273 with 
INTEGRAL* , XMM-Newton and RXTE 
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Abstract. INTEGRAL has observed the bright quasar 3C 273 on 3 epochs in January 2003 as one of the first observations of 
the open programme. The observation on January 5 was simultaneous with RXTE and XMM-Newton observations. We present 
here a first analysis of the continuum emission as observed by these 3 satellites in the band from ~ 3 keV to ~ 500 keV. The 
continuum spectral energy distribution of 3C 273 was observed to be weak and steep in the high energies during this campaign. 
We present the actual status of the cross calibrations between the instruments on the three platforms using the calibrations 
available in June 2003. 
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1. Introduction 



INTEGRAL ( Winkler et al. 2003) observed the bright quasar 
3C 273 (see Courvoisier 1998 for a review of the properties of 
the quasar) at 3 epochs in January 2003. This is the beginning 
of a programme meant to monitor the high energy emission of 
the object in order to measure the different components con- 
tributing to the emission above 1 keV and their respective vari- 
ability. Some observations of the INTEGRAL programme are 
conducted in coordination with the XMM-Newton and RXTE 
satellites in order to provide cross calibrations between the in- 
struments on board these three platforms. 

We report here a first analysis of the INTEGRAL data to- 
gether with the results of the XMM-Newton and RXTE observa- 
tions that were conducted simultaneously with the INTEGRAL 



measurements in January 2003. The INTEGRAL data span 12 
days while the XMM-Newton and RXTE data were obtained on 
a single date during the first INTEGRAL observations. 
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2. INTEGRAL observations 

INTEGRAL observed 3C 273 during revolution 28 starting 
January 5, 2003 for 1.2 x 10^ s, during revolution 30 starting 
January 11, 2003 for 10^ s and during revolution 32 starting 
January 17, 2003 for 1.1 x 10^ s. The observations were per- 
formed using the 25 points square dithering pattern that is ex- 
pected to give the best results for the SPI analysis except for 
23 000 s during which the satellite pointed stably on the source. 
Only the revolution 28 INTEGRAL observation was strictly si- 
multaneous with the RXTE and XMM-Newton observations. 
There is as of now no evidence of significant 3C 273 flux varia- 
tions observed by the INTEGRAL instruments between the rev- 
olutions 28, 30 and 32. We therefore combine here all available 
data. Future analysis will be done to quantify possible flux vari- 
ability. 
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2.1. IBIS/ISGRI results 

We present in Fig. 1 a significance image of the IBIS/ISGRI 
(Lebrunet al. 2003 1 25-40 keV data accumulated over 1.62 x 
10^ s. We used version 2.0 of ISDC's ( Courvoisier et al. 2003 1 
Offline Science Analysis (OSA) software. The algorithms used 
in the analysis are described in Goldwurm et al. 2003 3C 273 
and the bright Seyfert 2 galaxy NGC 4388 are clearly detected 
at a significance level of 10.5 for 3C 273 and 9.3 for NGC 4388. 
No source is detected in the surrounding of 3C 273. Other max- 
ima in the significance image will be analysed to assert their na- 
ture as either noise or possible sources below the present level 
of detection. The 3C 273 countrate in the 25 to 40keV band 
is 2.66 ± 0.12 counts s^^ (here and in the following we use 
1 <T confidence level when quoting uncertainties). The flux of 
3C 273 at 50keV is about 7mCrab. 

We have used the image shown in Fig. 1 to assess 
the reality of the existence of the source GRS1227H-025 
(Jourdainet al. 1992 1 which was claimed to be detected dur- 
ing a 1990 SIGMA observation. We show in the insert of Fig. 1 
a flux image of the region near 3C 273 with the GRS 1227H-025 
possible localisation marked. GRS1227H-025 was found in the 
1990 SIGMA data to be much brighter than 3C 273 which was 
not seen during that observation. [Leach and McHardy 1996| 
have searched for a possible counterpart to GRS1227H-025 in 
the ROSAT data and concluded that there is no likely candi- 
date and that the source, if real, must be transient or highly 
absorbed. Our data show no source at the position of GRS 
1227H-025 (Fig. 1). At the present level of the analysis we 
can conclude that a conservative upper limit to the flux of 
GRS1227H-025 is 1/2 that of 3C 273 in the 25-40 keV band. 
This estimate is based on our understanding of the system- 
atic uncertainties rather than on statistics. We expect that more 
stringent limits will be obtained in the future. The other source 
in the field of view that is seen in Fig. 1 is the bright Seyfert 



Fig, 1, IBIS/ISGRI significance image in the 25-40 keV energy 
band of the field of 3C 273 and NGC 4388. The insert shows 
a 25-40 keV flux image of the area surrounding 3C 273 and 
where GRS I227H-025 would be expected. The 3C 273 count 
rate is about 2.7 counts s~^, that of GRS 1227H-025 is less than 
0.4 counts s^^. 



2 galaxy NGC 4388. This source is 10.6° away and does not 
contribute therefore to the fluxes given by the instruments used 
in this study. 

We extracted the ISGRI spectrum of 3C 273 accumulated 
over 1.62 x 10^ s using the data from the revolutions 28, 30 and 
32. We used an intermediate version of the OSA software and 
the matrices describing the instrument as available at the ISDC 
in June 2003. We have ignored the channels below 25 keV and 
above lOOkeV. The resulting count spectrum can be well fitted 
with a single power law function between 25 and lOOkeV. The 
photon index is 1.95 ± 0.2. 

2.2. SPI results 

The analysis of the SPI (IVedrenne et al. 2003> data is based on 
74 dithering pointings with a total exposure time of 147 ksec. 
From the total of 83 dithering pointings which were taken dur- 
ing the 3C 273 observation, 9 had to be excluded from the SPI 
analysis as they either were affected by strong solar activity or 
had been influenced by the radiation belts. As the SPI data are 
background dominated, a careful background substraction is 
essential in order to get reasonable results, especially for weak 
sources. A time dependent background model has been applied 
to the data, based on the saturated events seen by the detec- 
tor The image reconstruction used for the analysis presented 
here is done by using the Iterative Removal Of Sources (IROS) 
method ( Hammersl ey et al. 1992^ which is im plemented in the 
SPIROS software iSkinner & ConneU 20031 . To get precise 
flux values, the source positions of the three brightest sources 
in the field (3C 273, 3C 279, and NGC 4388) have been fixed to 
their catalogue values, and fluxes have been extracted. No other 
source with a significance larger than 3a has been detected in 
the SPI data. For extraction of a rough spectrum, five energy 
bins (with boundaries at 20, 40, 100, 200, 500, 1000 keV) have 
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Fig. 2. OMC light curve. IJD is the fractional number of days 
since January 1, 2000 at UT 
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been applied to the data. 3C 273 is detectable in the SPI data up 
to at least 200 keV. The flux value in the 200 - 500 kcVband is 
/ = (5.5 ± 4.9) X lO""' photons cm~^ sec~^. Simulations of 
weak sources have also shown that the analysis software tends 
to overestimate their flux at high (> 200 keV) energies, and 
thus a detection up to 500 keV cannot yet be confirmed by the 
SPI data. The SPI spectrum is consistent with a single power 
law with photon index F — 1.66 ± 0.28 and the flux averaged 
over the whole observation is ~ 10.0 ± 2.5 mCrab. 

2.3. JEM-X results 

We analysed all the data in which 3C 273 was on-axis, this 
amounts to 32000 s of JEM-X 2 JLund et al 2003 1 data (the 
time during which INTEGRAL was in staring mode plus those 
pointings where the source was on-axis). The analysis was 
performed by fixing the position of the source, extracting the 
counts spectra for each pointing individually and adding the 
results weighted by the exposure time of the individual point- 
ings. 

The spectral extraction is based on an algorithm which 
is similar to photon tagging ("F enimore 1986> . Each detected 
photon is backprojected through the mask and its contribu- 
tion to the source and background flux estimates for the par- 
ticular source is computed. The spectrum is built up by the 
net contribution from all detected photons. The implemen- 
tation of the algorithm in the JEM-X case is described in 
IWestergaard et al. 2003| 

When performing a single power law spectral analysis of 
the extracted spectrum we obtain a slope F = 1.6 ± 0.23. The 
3-lOkeV flux is (4.0 ± 0.2) x 10"" ergscm^^ g-i. 

2.4. OMC light curve 

The Optical Monitoring Camera (OMC, 
IMas-Hesse et al. 2003> provides a simultaneous set of 
observations of several sources in the INTEGRAL field of 
view. We considered all exposures of 100 s of 3C 273 and ex- 
tracted square boxes (187 x 187arcsec^) centered on 3C 273. 
The images of these boxes have been combined in a single 
averaged image per science window and their analysis was 
then performed by "square photometry" using the IRAF image 
package IMSTATISTICS, first on the whole box and then on 
small boxes centered on sources close to 3C 273 in order to 
remove their contribution to the flux of the quasar and to the 
estimate of the background. Indeed as the pixel size of OMC 
CCD is ~ 17 arcsec and severly undersamples the point spread 
function (PSF), a circular aperture photometry or PSF fit are 
inapropriate. 

The resulting light curve (Fig. 2) shows an average flux 
level of my = 12.72 ± 0.02 during revolution 28, my — 
12.79 ± 0.03 during revolution 30 and my = 12.61 ± 0.03 
during revolution 32. Figure 2 shows some level of variation 
(~ 0.3 magnitude) during all revolutions. A careful examina- 
tion of the data gives no indication that instrumental effects 
could affect these measurements. In particular, these variations 
are still present when the light from the whole window is con- 



sidered. Furthermore, the 2 other sources in the window are too 
weak (V = 13.2 and V = 14.9) to account for the observed vari- 
ations. No source on the outside of the box is expected to influ- 
ence the total flux measured, the background was not observed 
to vary significantly. Further work is in progress to understand 
these variations. 

3. XMM-Newton Observations 

XMM-Newton observed 3C 273 January 5, 2003 from 14:22 to 
16:45 and from 17:29 to 18:53 (UT). 

From the EPIC MOS/PN cross-calibration' it is known that 
for the Small Window (SW) mode the spectral shape obtained 
with MOS2 and PN are consistent with one another. The MOS 1 
camera was in several different modes during the observations. 
The MOS2 camera was in the Small Window mode throughout 
the observation, its data is therefore retained for the remainder 
of this analysis. 

The MOS2 Small Window observations show pile-up in the 
region of the source center. Therefore the central region of the 
image is excluded from the analysis. To examine the source 
area containing pile-up we excluded a circular region of the 
source center and decreased the radii of this regions until pile- 
up effects become perceptible. This is verified by the compari- 
son of the pattern distribution and pattern fraction of valid pat- 
terns of type single, double, triple, and quadruple with model 
curves using the SAS V.5.4.1 task EPATPLOT V.1.1.8. 

The PN observations do not show pile-up effects. The small 
discrepancy seen in the pattern distribution below 600 eV is re- 
lated to a PN calibration issue and not related to pile-up. 

For the analysis presented in this paper circular source re- 
gions are used with size of 0^0 arcsec for PN SW and 10- 
40 arcsec for MOS2 S W. The spectra are generated out of sin- 
gle events (pattern 0) only. The PN background was selected 
from two boxes near the border of the small window chip range. 
For M0S2 it is impossible to define a background region on 
CCDl because the PSF fills completely the small window area. 
The background is determined from two circular regions of 
CCD3/6 where no background source is identifiable. 

We performed a spectral fit above 3 keV, thus avoiding the 
complex soft X-ray excess region, for the first period of ob- 
servation. We used a simple power law and obtain a slope 
r = 1.74 ± 0.03 (PN) and 1.73 ± 0.08 (MOS2). The 3- 
lOkeV fluxes are (6.8 ± 0.7) x 10"" ergs cm^^ g-i (PN) and 
(7.0 ± 1.8) X 10-"ergscm-2s-i (MOS2). The flux uncer- 
tainties take the uncertainty on the slope and the normalisation 
into account. These results are entirely consistent with those 
obtained during the second period of observation. 

There is no evidence of an Fe Kq line. Upper limits for the 
equivalent width of a neutral or ionised Fe line are < 52 eV 
for a narrow (cr = 0.01 keV) and < 90 eV for a broad (cr = 
0.5keV)line. 

4. RXTE Observations 

RXTE observed 3C 273 January 5, 2003 between 14:22 and 
19:00 (UT). 

' XMM-SOC-CAL-TN-0018, 2003-04-04 
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Fig. 3. The spectral energy distributions measured by the 3 high energy instruments of INTEGRAL, by XMM-Newton and 
RXTE. 



For RXTE only PCUO and PCU2 were enabled. Because the 
count rate of the quasar observation was very low only the top 
layers were taken into account to maximise the signal to noise 
ratio. Since May 2000 the top layer of PCUO is damaged, so 
only the top layer of PCU2 was used for spectral analysis. 

The PCA Standard Mode 2 data were reduced us- 
ing FTOOLS V5.2. The response file was generated using 
PCARMF V8.0. The background estimation was done with 
PCABACKEST V3.0 using the L7-240 background model rec- 
ommended for faint sources by the NASA/GSFC RXTE guest 
observer facility. 

Both HEXTE clusters were combined to get higher signal 
to noise ratio. Nevertheless the count rate is very weak. The 
poor statistics causes PCA/HEXTE cross-calibration problems. 
Fitting the PCA and HEXTE (20-190 keV) data simultane- 
ously the fit is completely dominated by the PCA, the result- 
ing r(PCA/HEXTE) is identical to r(PCA). We show here the 
HEXTE data over the complete energy range, aware of its low 



significance, to show how it compares with the other instru- 
ments. 

The results of the PCA spectral analysis give F — 1.73 ± 
0.02anda3-10keV flux of (9.2±0.5) x 10"" ergscm^^ g-i. 
There is no evidence for an Fe line. 



5. Discussion 

We present in Fig. 3 the spectral energy distributions obtained 
by the high energy instruments of INTEGRAL and by XMM- 

Table 1. Intercalibration factors deduced from a fit to the 
3C 273 high-energy data with a power law. 



XMM-PN 


XMM-MOS2 


RXTE-PCA 


RXTE-HEXTE 


1.0 (fix) 


1.024 ±0.025 


1.34 ±0.02 


0.93 ± 0.06 


JEM-X 


IBIS-ISGRI 


SPI 




0.61 ±0.05 


1.76 ± 0.2 


0.95 ±0.3 
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Fig. 4. The spectral energy distributions of figure 3 normalised to the PN camera of XMM-Newton flux and compared with the 
historic average (continuous line) and observed range of variations as given in lTiirler et al. 19991 



Newton and by RXTE. The photon distributions as well as a 
V fu representation are given. The data are shown without tak- 
ing into account the factors that correct the normalisations ob- 
tained with the different instruments. It is clear from Fig. 3 and 
Table 1 that the measurement of the detailed shape of the emis- 
sion over large spectral ranges spanning more than one instru- 
ment as it is needed to establish the presence of reflection com- 
ponents and cut-offs hinges on the intercalibrations between 
the instruments not only on different satellites but also on the 
same platform. 

The average V magnitude of 3C 273 quoted in 
ICourvoisier 1 998 is 12.9 mag. The magnitudes measured here 
are between 12.61 mag and 12.79 mag showing that the object 
was slightly brighter in the optical domain during this cam- 
paign than it usually is. 

The high energy continuum spectral energy distribution 
measured here (Figure 4) using cross calibration factors of the 
instruments given in Table 1 (i.e. aligning all instruments to the 



PN Camera of XMM-Newton) is well represented by a power 
law with a photon index of 1.73 ± 0.015 and normalisation of 
(2.24±0.05) X 10-2 photons cm'^ s"! keV^i at 1 keV. The re- 
duced chi-square is 0.94 for 292 degrees of freedom. This mea- 
surement differs from the spectrum of lLichti et al. 19951 which 
has a similar normalisation at 1 keV ((2.33 ±0.04) x 10"^ pho- 
tons cm^^ s^^ keV^^) but is flatter, the slope being 1.6 ± 0.01. 

The continuum spectral energy distribution shown in Fig. 
4 is normalised to the PN camera of XMM-Newton flux value. 
It is slightly below the historic average of lTiirler et al. 1999l in 
the low energies. The slope deduced from the measurements 
presented here is steeper than that suggested by the historical 
average which is close to that of L ichti et al. 1995 The contin- 
uum spectral energy distribution between 10 and 100 keV is 
therefore close to the weakest observed to date (Fig. 4). 

We conclude that there is an important cross calibration 
effort to be performed in order to be able to confidently use 
multi-instrument data to model the emission of sources over 
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several decades of energy. Fig. 3 indicates that there is a sig- 
nificant mismatch in the present caHbrations of ISGRI and SPI 
(June 2003). The normalisation of 3C 273 deduced with SPI 
(Table 1 ) agrees better with that obtained from XMM-Newton 
than the normalisation obtained with ISGRI. The relative nor- 
malisations obtained with both RXTE instruments indicates that 
some cross calibrations issues remain to be solved. The exten- 
sive SPI ground calibration campaign (lAttie et al. 2003t will 
prove very important in settling the issue of cross-caUbration 
in the near future. 

3C 273 has been weak in the high energy domain in January 
2003 and the spectral slope is steeper than in previous cam- 
paigns. There is no evidence for components other than a single 
power-law above 3 keV in the data presented here. 
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